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Abstract 



We investigate the critical dynamics of the time-dependent Landau-Ginzburg 
model with non conserved n-component order parameter (Model A) in the 
presence of long-range correlated quenched impurities. We use a special kind 
of long-range correlations, g{\r — r'\) ~ |r — r'\~°', a = 4 — 5, previously 
introduced by Weinrib and Halperin. The parameter (5 ~ e, where e = A — d 
is the usual small parameter within the Wilson-Fisher renormalization group 
theory. 

Using a double expansion in e and 5 we calculate the critical exponent 
z up to second order on the small parameters. We show that the quenched 
impurities of this kind affect the critical dynamics already in first order of e 
and 5, leading to a relevant correction for the mean field value of the exponent 
z. 
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The static and dynamic scaling properties of isotropic systems with short-range (SR) 
correlated quenched impurities has been the subject of a large number of studies, and is now 
well understoodEhi. From those works, it is known that the critical behavior, described by the 
SR fixed point (f.p.) of the renormalization group (RG) recursion relations, is consistent with 
the Harris criterion!, which states under which conditions a random critical-temperature 
system undergoes a second order phase transition with the same critical properties as the 
corresponding pure system. 

For the static proDerties, the case of long-range (LR) correlations was considered by 
Weinrib and HalperinQ. They introduced a special kind of spatial Gaussian correlations 
between the quenched impurities falling off with the distance as a power law. Their model 
permits to find a new f.p. of the RG transformation corresponding to a second order phase 
transition with specific critical exponents due to the presence of this kind of impurities. 

When disorder with LR power-low correlations is taken into account, the Harris criterion 
is modified. If a is the exponent of the LR correlations (see below), d is the dimensionality 
of the system, and v is the correlation- length exponent ^ ~ |T — Td^*^, the modified Harris 
criterion! states that the LR disorder is irrelevant if az/ — 2 > for a < ci, while dv ~ 2 = 
—a > for a > d (the normal Harris criterion). 

As a further step in that direction, in this work we analyze, in the framework of the 
e-expansion, the critical dynamics of the time-dependent Landau- Ginzburg (TDLG) model 
for systems with n-component , non-conserved order parameter (Model A)ll3, in the case 
when the LR correlations between impurities are of the form introduced by Weinrib and 
Halperin. We calculate the dynamic critical exponent z up to second order in e = du — d 
{d is the dimensionality of the space and (i„ = 4 is the upper borderline dimensionality for 
this model), generalizing the result of Grinstein, Ma and MazenkoB, valid for the case of SR 
correlation between impurities. 

The model we consider is 

dipa{x,t) dF 

where ipa{x,t),a = 1, n is a n-component order parameter, F is a bare kinetic coefficient, 
Tja is a Gaussian white noise 

<rj^{x,t) >av=0, (2) 



< ria{x,t)rif3{x\t') >av= 2T5{x - x')S{t - t')5a,f3 



(3) 



and 



d'^x 



rip^x, t) + (yij{x, t)f + -u(?/'2(f, t)f + (^(f)?/'2(x, t) 



(4) 



is the Landau-Ginzburg free energy functional including a static random term V9(x), which 
describes the quenched impurities. 

This term obeys a Gaussian distribution 



< V9(x) > 



av 0) 



(5) 
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< Lp{x)(f{x') >av= g{\x - X'\), (6) 

where the function g{r), r = \x — x'\ , is a. hnear superposition of inverse power law correla- 
tions decaying with the distance r at different rates g{r) ~ Si^i Air'""*. It is known! that 
of the different components of g{r), the most relevant for the scaling analysis is that with 
the slowest decay rate. 

The Fourier transform of g{r) for small k iM 

g{k) = V + Wk^-^. (7) 

The introduction of a second small parameter 5 = 4 — a, in addition to e, 5 = 0(e), 
and the double expansion in e and 5 permits to studv the nontrivial LR critical behavior 
described by a new f.p. of the RG recursion relational. This important result is due to the 
similar scaling behavior when 5 ~ 0{t) of the two presumably different (short- and long- 
range) behaviors implied in (|^. Actually, the inequality a > d leads in the long- wavelength 
limit {k — > 0) to the scaling irrelevance of the parameter w and of any effect of the LR 
correlations. When a < d, the second term in (^ is dominant for small k, but leading to 
an unstable critical behavioiS. Both terms have the same scaling relevance for the critical 
behavior only if a ~ , which will be the case considered herel. 

The calculation of the critical exponent z closely follows ref.@. Using the TDLG equations, 
(|^-^) we calculate the renormalized response function from the Dyson equation, looking for 
self-energy o^-dependent contributions of the form {—iuj) ln(— ■ju;)llHl3. A trivial calculation 
leads to the following expression for the dynamic exponent z with a correction already in 
first order in e and S 

z = 2 + {v* + w*), (8) 

where v* and w* are the f.p. values corresponding to the LR f.p. for the parameters v and 
w, derived from the recursion relations to order 0(e)li 

^ 2(n + 2)e- (n + 8)(5)', (9) 



2(5-e)(5r2 + 4)2 



= 2(^-e)(5n + 4)2 ^^"^ + " ^^"^ + ^^'^ ^^^"^ " " ^""'^ • 
For completeness we give here the f.p. value of u, which will be used later 

(5n + 4) ^ ^ 

After substituting the f.p. values v* and w* in Eq.(|^) we get 

^ = 2+ „.- ^ ((n + 8)(5-2(n + 2)e), (12) 
2[bn + 4j 

valid up to order 0(e, S). 

Notice that an e-correction for z is also present in the SR ca.s(^ 
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The calculation of the second-order correction in e for the dynamic exponent z is analo- 
gous to that ini. We again perform the analysis as before by calculating the {—ioj) ln(— ia;) 
and {—ioS) ln^(— icu) contributions up to order e^. To this order the response function can be 
written as 



-iuj\ , ( -ioj 
- l + ln — 



^ + 2)m^ n iv ^wY ^ 

2 8^ ^ 3 2^ ^ 4 



+ ln' 



v + w\ 5, ,2 n + 2 . . 1 , ^' 
e + -{v + w) — u[v + w) — -w[e — o) 



-u{v + w 

(13) 



Note that the contribution —w{e — 5)/8 coming from the diagram proportional to f + tt; 
ensures that the relation 



1 /f + ^ 

2 V~2" 



V + w 



e + -{v + wY - 



n+2 . , 1 , 

— - — u[v + w) — -w{e — o) 



(14) 



is consistently satisfied by the 0(e) f.p. values of the parameters u^v and w (Eqs.p|-|Tl|). 
From the response function we obtain the expression for z up to order 0(e^, e5, 5"^) 



z = 2-ri + {v + w) + -in + 2)u^ hi- + -{v + wf - ^ ' 



3 

4' 



4 3 
3 ■ 2* 



u{y + w), 



where the anomalous dimension r] ialll 

1 



V 



16(5n + 4)^ 



8n{n + 2)e' + {n + 4)^5(2e - 36) 



(15) 



(16) 



From (|T5D it is clear that in order to calculate the dynamic exponent z to order e^, one 
needs the expression up to this order only for the f.p. value of the term v + w. For the last 
three terms, which are quadratic, it is enough to use the f.p. values of the parameters given 
by (1^0). 

We calculate the f.p. values up to order by resolving the following recursion relations 
for M, V and w written to that order 



u' = 6^-2" \u- Kd\nb[l + - \nb 



u — bu{v + w) 



+Kl\nH 



1 3 

-(n^ + Qn + 20)u^ - -{n + 8)u'^{v + w) + 9u{v + wf 



+ Kl\nh{l + \iih) 



■{5n + 22)u^ - 6(ra + ^)u^{v + w) + 2lu{v + wf 



(17) 



v' = b'-^"^ \ v - Kdlnb{l + - \nb) \ {n + 2)uv - 2v{v + vS) - 2{v + wY\ - 3vw{e - b)K^ In^ b 



-2nY[e - 5)Ki In^ b+Kj In^ b 
+7^4^ In 6(1 + In 6) 
— 3{n + 2)uv{v + w) 



2{v + wf + 3v{v + wf + -{n + 2fu^v - 3{n + 2)uv{v + w) 
3 

{v + wf + 3v{v + wf + -{n + 2)u% - 3{n + 2)u{v + wf 

2 
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w' = 6^-2" S^w - Kdlnb (l + ^ In [{n + 2)uw -2w{v + w)] - w^(e - 6)K^ In^ b 



+Kj In' b 



3w{v + wY + -(n + 2Yu^w — ?>{n + 2)uw{y + w) 



ir42ln6(l + ln6) 



?>w{y + w)^ + -(n + 2,)u^w — 3(n + 2)uw{v + w) 



(19) 



In (17-19), 6>1 is a rescaling factor, and all integrals over momenta were taken within 
the shell b^^ < k < 1. The constant = 2^('^^^)7r~'^/2/r(|) will be later absorbed in the 
definition of the f.p. values u*, v* and w*, in correspondence with (|9|-[TTD. Notice also that 
for deriving the above recursion relations we ignored the mass-renormalization graphsEBIii. 
Putting w = in (p!7|-p!9|) one immediately obtains the recursion relations for the SR caseSi. 

After some algebra, including correctly all terms up to order e^, one can see that all the 
In 6 contributions in the recursion relations cancel when one calculates the f.p. values of 
the parameters, leading to the following expression for v* + w* (here we write only its value 
because of the large and complicated f.p. expressions for the other parameters) 



v* + w*= \ {in + 8)5 -2{n + 2)e) + . . . .3 
z(5n + 4) lo(5n + 4)'^ 



9(15^3 + 76n2 + 16n + 64)5^ 



+ 2(-95n^ - 108^2 + 240n - 64)e5 - 32(n + 2f{bn + 2)e^ 

Substituting this last expression in (1T5|), we finally obtain up to order 0(e^, e5, 5^ 

{n + 2) 



(20) 



2 + 



1 



2(5n + 4) 
4(10^2 + 19n + 4)e2 



((ri + 8)5-2(n + 2)e) + 



4(5n + 4)3 



3(n-4)(5n-8)(5(5- -e) 



+ 



3(n + 2) 
(5n + 4)2 



(5) 









(21) 



which is valid for n > 1. The case n = 1 should be treated separately because of the 
degeneracy of the recursion relationdli'i. According toi, the inclusion of the term w 7^ 
breaks the degeneracy, and in order 0(e) the LR disorder f.p. found before is a solution of 
the recursion relation for n = 1. It is however unstable for 5 ~ Oie). Performing thus an 
expansion in e^l\ and assuming 5 = 0(eV2), one obtains for the f.p. values! 



u 



0(e), 



(22) 



which leads to results consistent with the extended Harris criterion. 

With the aid of (|2^), we get from (|T5]) the dynamic exponent for the case n = 1, 

z = 2^'- 



0(e) 



(23) 



Equations (pTD and ( p3D complete our study of the critical dynamics of the TDLG equa- 
tions in the presence of LR correlated quenched impurities. A further investigation could 
be the study of the regions of stability of the LR critical behavior taking into account the 
e^corrections. This, however, will not affect the qualitative behavior of the phase transition, 
and a similar picture to that described inll, characterized by oscillatory corrections to the 
scaling, will appear. 
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In summary, in this work we investigated the dynamic scahng of the TDLG equations 
with non-conserved order parameter (Model A) in the presence of quenched impurities with 
LR spatial correlations. Using a double expansion in two small parameters e and 5 ^ 0{e), 
we calculated the critical exponent z up to order O(e^). We showed that the quenched 
impurities of this kind affect the critical dynamics already in first order in e and 5, which 
is a common characteristic for systems with quenched impurities, independently of the kind 
of correlations (short- or long-range). 

This work is supported by DGICYT project PB94-0388, NATO grant CRG95055 and 
partially by Contract F-608 from the Bulgarian Scientific Foundation. 
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